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' Fractal structures are observed in the universe in two very different ways. 

' Firstly, in the gas forming the cold interstellar medium in scales from 10~^pc 

o 



o 



till lOOpc. Secondly, the galaxy distribution has been observed to be fractal 
00 '. in scales up to hundreds of Mpc. We give here a short review of the statisti- 

cal mechanical (and field theoretical) approach developed by us for the cold 
interstellar medium (ISM) and large structure of the universe. We consider 
a non-relativistic self-gravitating gas in thermal equilibrium at temperature 
T inside a volume V. The statistical mechanics of such system has special 
^ , features and, as is known, the thermodynamical limit does not exist in its 

customary form. Moreover, the treatments through microcanonical, canoni- 
cal and grand canonical ensembles yield different results. We present here for 
^ ' the first time the equation of state for the self-gravitating gas in the canon- 

■ ical ensemble. We find that it has the form p = [NT/V]f{ri), where p is 

the pressure, is the number of particles and ri = ^y/sy- '^^^ N ^ oo 
and V ^ oo limit exists keeping r] fixed. We compute the function /(rj) 
using Monte Carlo simulations and for small r], analytically. We compute 
the thermodynamic quantities of the system as free energy, entropy, chemical 
potential, specific heat, compressibility and speed of sound. We reproduce 
the well-known gravitational phase transition associated to the Jeans' insta- 
bility. Namely, a gaseous phase for rj < rjc and a condensed phase for rj > rjc- 
Moreover, we derive the precise behaviour of the physical quantities near the 
transition. In particular, the pressure vanishes as p ~ {rjc — rf)^ with B ~ 0.2 
and ric ~ 1.6 and the energy fluctuations diverge as ~ (r/c — v)^^^- The speed 
of sound decreases monotonically with r/ and approaches the value \JT at 
the transition. 
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I. STATISTICAL MECHANICS OF THE SELF-GRAVITATING GAS 



Physical systems at thermal equilibrium are usually homogeneous. This is the case 
for gases with short range intermolecular forces (and in absence of external fields). When 
long range interactions as the gravitational force are present even the ground state is often 
inhomogeneous. In this case, each element of the substance is acted on by very strong 
forces due to distant particles of the gas. Hence regions near to and far from the boundary 
of the volume occupied by the gas will be in very different conditions and as a result the 
homogeneity of the gas is destroyed |^. This basic inhomogeneity suggests that fractal 
structures can arise in a self-interacting gravitational gas 

Let us review very briefly our recent work [0-^ on the statistical properties of a self- 
interacting gravitational gas in thermal equilibrium. We discussed two relevant astrophysical 
applications of such gas: the cold interstellar medium (ISM) and the galaxy distributions. 

In the grand canonical ensemble, we showed that the self-gravitating gas is exactly equiv- 
alent to a field theory of a single scalar field 0(a;) with exponential self-interaction. We 
analyzed this field theory perturbatively and non-perturbatively through the renormaliza- 
tion group approach. We showed scaling behaviour (critical) for a continuous range of the 
temperature and of the other physical parameters. We derive in this framework the scaling 
relation 

M{R) ~ R'^» 

for the mass on a region of size R, and 

Av ~ i?'? 

for the velocity dispersion where q = ^{dn — 1). For the density-density correlations we find 
a power-law behaviour for large distances 

The fractal dimension dn turns to be related with the critical exponent u of the correlation 
length by 

dH= l/u . 

Mean field theory yields for the scaling exponents u = 1/2, dn = 2 and q = 1/2. Such values 
are compatible with the present ISM observational data: 1.4 < dn 0.3 < q < 0.6 . 

We developed in ref. [Q a field theoretical approach to the galaxy distribution. We 
consider a gas of self-gravitating masses on the Friedman- Robertson- Walker background, in 
quasi-thermal equilibrium. We derive the galaxy correlations using renormalization group 
methods. We find that the connected A^-points density correlator C{ri,r2, ■ ■ ■ , r/v) scales as 

N{D-3) 
'1 ) 

when ri >> Tj, 2 < i < N. There are no free parameters in this theory. 

Our study of the statistical mechanics of a self-gravitating system indicates that gravity 
provides a dynamical mechanism to produce fractal structure. 

This paper is organized as follows. In section II we summarize the main properties of the 
ISM, in section III we review the relevant aspects of the large scale structure of the universe, 
in sec. IV we develop the statistical mechanics of the self-gravitating gas in the canonical 
ensemble. Discussion and remarks are presented in section V. 
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II. THE INTERSTELLAR MEDIUM 



The interstellar medium (ISM) is a gas essentially formed by atomic (HI) and molecular 
{H2) hydrogen, distributed in cold (T ~ 5 — 507^) clouds, in a very inhomogeneous and 
fragmented structure. These clouds are confined in the galactic plane and in particular 
along the spiral arms. They are distributed in a hierarchy of structures, of observed masses 
from 10~^ Mq to IO^Mq. The morphology and kinematics of these structures are traced 
by radio astronomical observations of the HI hyperfine line at the wavelength of 21cm, 
and of the rotational lines of the CO molecule (the fundamental line being at 2.6mm in 
wavelength), and many other less abundant molecules. Structures have been measured 
directly in emission from O.Olpc to lOOpc, and there is some evidence in VLBI (very long 
based interferometry) HI absorption of structures as low as 10~^ pc = 20 AU (3 10^^ cm). 
The mean density of structures is roughly inversely proportional to their sizes, and vary 
between 10 and 10^ atoms /cm? (significantly above the mean density of the ISM which is 
about 0.1 atoms/cm^ or 1.6 10"^^ gjcm? ). Observations of the ISM revealed remarkable 
relations between the mass, the radius and velocity dispersion of the various regions, as first 
noticed by Larson 0, and since then confirmed by many other independent observations 
(see for example ref. 0). From a compilation of well established samples of data for many 
different types of molecular clouds of maximum linear dimension (size) i?, total mass M and 
internal velocity dispersion Af in each region: 

M{R) ~ R^^ , At; ~ , (2.1) 

over a large range of cloud sizes, with 10~^ — 10^^ pc < R < 100 pc, 

1.4 <dH<2,0.3<q< 0.6 . (2.2) 

These scaling relations indicate a hierarchical structure for the molecular clouds which is 
independent of the scale over the above cited range; above 100 pc in size, corresponding to 
giant molecular clouds, larger structures will be destroyed by galactic shear. 

These relations appear to be universal, the exponents (Ih, q are almost constant over all 
scales of the Galaxy, and whatever be the observed molecule or element. These properties 
of interstellar cold gas are supported first at all from observations (and for many different 
tracers of cloud structures: dark globules using ^^CO, since the more abundant isotopic 
species ^^CO is highly optically thick, dark cloud cores using HCN or CS as density tracers, 
giant molecular clouds using ^^CO, HI to trace more diffuse gas, and even cold dust emission 
in the far- infrared) . Nearby molecular clouds are observed to be fragmented and self-similar 
in projection over a range of scales and densities of at least 10"^, and perhaps up to 10^. 

The physical origin as well as the interpretation of the scaling relations ( |2.1| ) are not 
theoretically understood. The theoretical derivation of these relations has been the subject 
of many proposals and controversial discussions. It is not our aim here to account for all 
the proposed models of the ISM and we refer the reader to refs. [0 for a review. 

The physics of the ISM is complex, especially when we consider the violent perturbations 
brought by star formation. Energy is then poured into the ISM either mechanically through 
supernovae explosions, stellar winds, bipolar gas flows, etc.. or radiatively through star 
light, heating or ionising the medium, directly or through heated dust. Relative velocities 
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between the various fragments of the ISM exceed their internal thermal speeds, shock fronts 
develop and are highly dissipative; radiative cooling is very efficient, so that globally the ISM 
might be considered isothermal on large-scales. Whatever the diversity of the processes, the 
universality of the scaling relations suggests a common mechanism underlying the physics. 

We propose that self-gravity is the main force at the origin of the structures, that can be 
perturbed locally by heating sources. Observations are compatible with virialised structures 
at all scales. Moreover, it has been suggested that the molecular clouds ensemble is in 
isothermal equilibrium with the cosmic background radiation at T ~ ?)K in the outer parts 
of galaxies, devoid of any star and heating sources [^]. This colder isothermal medium might 
represent the ideal frame to understand the role of self-gravity in shaping the hierarchical 
structures. Our aim is to show that the scaling laws obtained are then quite stable to 
perturbations. 

Till now, no theoretical derivation of the scaling laws eq.( p.l|) has been provided in which 
the values of the exponents are obtained from the theory (and not just taken from outside 
or as a starting input or hypothesis). 

The aim of our work H-pl is to develop a theory of the cold ISM. A first step in this goal 



was to provide a theoretical derivation of the scaling laws eq. (|2.1| ), in which the values of the 



exponents g are obtained from the theory For this purpose, we implemented for 

the ISM the powerful tool of field theory and the Wilson's approach to critical phenomena 



III. GALAXY DISTRIBUTIONS 



One obvious feature of galaxy and cluster distributions in the sky is their hierarchical 
property: galaxies gather in groups, that are embedded in clusters, then in superclusters, 
and so on [^,0 . 

The knowledge of the galaxy and cluster correlations allows a more precise characteri- 
zation of their distributions. Unfortunately, the most widely spread two point correlation 
function ^(r) in galaxy distributions studies, is based on the assumption that the Universe 
reaches homogeneity on a scale smaller than the sample size, ^(r) is defined as 



< n[ri).n[ri + r) > 

< n >2 
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where n(r) is the number density of galaxies, and < ... > is the volume average (over d^Ti). 
The length tq is defined by ^(^o) = 1. The function ^(r) has a power-law behaviour of slope 
—7 for r < ro, then it turns down to zero rather quickly at the statitistical limit of the 
sample. This rapid fall leads to an over-estimate of the small-scale 7. One finds the slope 
7, the same for galaxies and clusters, of ~ 1.7 (e.g. 



It has been shown in refs. and pO| that the homogeneity hypothesis could perturb 
significantly the results. 



Pietronero |]32| introduces the conditional density 



r(r) 



< n{ri).n{ri + r) > 



<n> 
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which is the average density around an occupied point. For a fractal medium, where the 
mass depends on the size as 

M(r) oc 

D being the fractal (Haussdorf) dimension, the conditional density behaves as 

r(r) oc r^~^ 

This is exactly the statistical analysis used for the interstellar clouds, since the ISM as- 
tronomers have not adopted from the start any large-scale homogeneity assumption (cf. 

m)- 

The fact that for a fractal the correlation ^(r) can be highly misleading is readily seen 
since 



<n> 

and for a fractal structure the average density of the sample < n > is a decreasing function 
of the sample length scale. In the general use of ^(r), < n > is taken for a constant, and we 
can see that 

D = 3-7 . 

If for very small scales, both ^(r) and r(r) have the same power-law behaviour, with the 
same slope —7, then the slope appears to steepen for ^(r) when approaching the length tq. 
This explains why with a correct statistical analysis |3^, the actual 7 ^ 1 — 1.5 is smaller 
than that obtained using ^{r). This also explains why the amplitude of ^(r) and tq increases 
with the sample size, and for clusters as well. 

This scale- invariance has suggested very early the idea of fractal models for the clustering 
hierachy of galaxies |p^p8| . Since then, many authors have shown that a fractal distribution 



indeed reproduces quite well the aspect of galaxy catalogs, for example by simulating a fractal 
and observing it, as with a telescope |]3^ , |37[ . 



There is some ambiguity in the definition of the two-point correlation function ^(r) above, 
since it depends on the assumed scale beyond which the universe is homogeneous; indeed it 
includes a normalisation by the average density of the universe, which, if the homogeneity 
scale is not reached, depends on the size of the galaxy sample. Once ^(r) is defined, one 
can always determine a length tq where ^{tq) =1 |2T|. For galaxies, the most frequently 
reported value is ro ~ bh~^ Mpc (where h = i^o/lOOkm s~^Mpc~^), but it has been shown 
to increase with the distance limits of galaxy catalogs p2[, ro is called 'correlation length' 
in the galaxy literature. [The notion of correlation length ,^0 is usually different in physics, 
where ^0 characterizes the exponential decay of correlations (~ e"*^/^"). For power decaying 
correlations, it is said that the correlation length is infinite]. 

The same problem occurs for the two-point correlation function of galaxy clusters; the 
corresponding ^ (r) has the same power law as galaxies, their length tq has been reported to 
be about tq ^ 25h^^ Mpc, and their correlation amplitude is therefore about 15 times higher 
than that of galaxies |33|. The latter is difficult to understand, unless there is a consider- 
able difference between galaxies belonging to clusters and field galaxies (or morphological 
segregation). The other obvious explanation is that the normalizing average density of the 
universe was then chosen lower. 



This statistical analysis of the galaxy catalogs has been criticized in refs. 32,26,20 



who stress the unconfortable dependence of ^{r) and of the length tq upon the finite size 
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of the catalogs, and on the a priori assumed value of the large-scale homogeneity cut- 
off. A way to circumvent these problems is to deal instead with the average density as a 
function of size. It has been shown that the galaxy distribution behaves as a pure self-similar 
fractal over scales up to ~ lOO/i^^ Mpc, the deepest scale to which the data are statistically 
robust p8[. This is more consistent with the observation of contrasted large-scale structures, 



such as superclusters, large voids or great walls of galaxies of ~ 200/i~ Mpc After a 
proper statistical analysis of all available catalogs (CfA, SSRS, IRAS, APM, LEDA, etc.. for 
galaxies, and Abell and AGO for clusters) Pietronero et al []32| , p8| state that the transition 
to homogeneity might not yet have been reached up to the deepest scales probed until now. 
At best, this point is quite controversial, and the large-scale homogeneity transition is not 
yet well known. 

Isotropy and homogeneity are expected at very large scales from the Cosmological Prin- 
ciple (e.g. 0]). However, this does not imply local or mid-scale homogeneity (e.g. p8|], ^8 



a fractal structure can be locally isotropic, but inhomogeneous. The main observational 
evidence in favor of the Cosmological Principle is the remarkable isotropy of the cosmic 
background radiation (e.g. [^), that provides information about the Universe at the mat- 
ter/radiation decoupling. There must therefore exist a transition between the small-scale 
fractality to large-scale homogeneity. This transition is certainly smooth, and might cor- 
respond to the transition from linear perturbations to the non-linear gravitational collapse 
of structures. The present catalogs do not yet see the transition since they do not look up 
sufficiently back in time. It can be noticed that some recent surveys begin to see a different 
power-law behavior at large scales (A ~ 200 — A00h~^ Mpc, e.g. ||27||). 

There are several approaches to understand non-linear clustering, and therefore the dis- 
tribution of galaxies, in an infinite gravitating system. Numerical simulations have been 
widely used, in the hope to trace back from the observations the initial mass spectrum 
of fiuctuations, and to test postulated cosmologies such as CDM and related variants (cf 
PP|). [That is numerically solving Newton's equations of motion of self-gravitating parti- 
cles]. This approach has not yet yielded definite results, especially since the physics of the 
multiple-phase universe is not well known. Also numerical limitations (restricted dynamical 
range due to the softening and limited volume) have often masked the expected self-similar 
behavior. 

We presented in a new approach based on field theory and the renormalisation group 
to understand the clustering behaviour of a self-gravitating expanding universe. We also 
consider the thermodynamics properties of the system, assuming quasi-equilibrium for the 
range of scales concerned with the non-linear regime and virialisation. Using statistical field 
theory, the renormalisation group and the finite-size scaling ideas, we determined the scaling 
behaviour. The small-scale fractal universe can be considered critical with large density 
fiuctuations developing at any scale. We derived the corresponding critical exponents which 
yielded the fractal dimension D. It is very close to those measured on galaxy catalogs 
through statistical methods based on the average density as function of size; these methods 
reveal in particular a fractal dimension D ?a 1.5 — 2 P3| , |S5[] . This fractal dimension is 
strikingly close to that observed for the interstellar medium or ISM (e.g. f^J^l) showed 
in ref. [^] that the theoretical framework based on self-gravity that we have developped for 
the ISM is also a dynamical mechanism leading to the fractal structure of the universe. 
This theory is powerfully predictive without any free parameter. It allowed to compute the 
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A^-points density correlations without any extra assumption 



IV. THERMODYNAMICS OF THE SELF-GRAVITATING GAS: THE 

CANONICAL ENSEMBLE 



We investigate in this section a gas formed by non-relativistic particles with mass 
m interacting only through Newtonian gravity and which are in thermal equilibrium at 
temperature T = We shall work in the canonical ensemble assuming the gas being on 
a cubic box of side L. 

The partition function of the system can be written as 

Z= [...fli'^f^e-^-^ (4.1) 



where 



t\ (2vr)= 



1=1 l<l<j<N W - 1j\ 



G is Newton's gravitational constant. 

Computing the integrals over the momenta pi, (1 < / < A^) 



J (2^ " 1^2^ J 



3/2 



yields 



Z 



m 



271(3 / Jo Jo 



r ... d^qi e^^"'^^^'<^^'^ l^^-'^^l (4.3) 

Jo Jo f-J^ 



We make now explicit the volume dependence introducing the new variables fi, 1 < I < N 

as 

qi = Lfi , fi = {xi,yi,zi) , 

0<xi,yi,zi<l . (4.4) 
The partition function takes then the form, 

Z = ^ ■■■ li d'n (4.5) 

\ 2tc J Jo Jo 

where we introduced the variable r], 

Gm'^N ,^ , 

^ = (4.6) 

and 
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M(fi,...,rjv) = 4 E 



Recall that 



Grri^N ,^ ^ , , , 

U = u{n,...,rN) (4.7) 



is the potential energy of the gas. 

The free energy takes then the form, 



F = -TlogZ = -SiVTlog I J^L I - T ^^{r]) (4.8) 



Here 

N 



$7v(r/) = log / . . . / n^^'"' 

"'0 Jo iJ[ 



•q u{fi,...,TN) 



(4.9) 



The derivative of the function ^n{v) will be computed by Monte Carlo simulations and, in 
the weak field limit rj << 1, it will be calculated analytically. 
We get for the pressure of the gas. 



[Here, V = L stands for the volume of the box]. We see from eg. ( |4.9| ) that ^Niv) increases 



with 7]. Therefore, the second term in eq. (|4.1CI|) is a negative correction to the perfect gas 



NT 

pressure 



The mean value of the potential energy U can be written from eq.( [4.7|) as 

< U >= -Tt] ^'j^{7]) (4.11) 
Combining eqs. ([4.1CI|) and ([4.11|) yields the virial theorem. 



or more explicitly 



where, 



pV _^ <U> 
iVT ~ 3NT 



NT 37V ^^'^ ^ ^ 



N 

^r)u{ri,...,rN) 



<l>'^(r/)=e-*-('') / .../ n^'nw(ri,...,r-W)e'^ 
Jo Jo iJ[ 

= -{N-l) e-*^^'') / . . . / n ^ ^ e''"^^!'-'"^^) (4.13) 

2 JO Jo ,1 Ti — r2 
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This formula indicates that ^'i^{f]) is of order for large A^. Monte Carlo simulations as 
well as analytic calculations for small rj show that this is indeed the case. In conclusion, we 
can write the equation of state of the self-gravitating g 

^ = m , (4.14) 

where the function 

fiv) = 1 - ^ ^Niv) , 

is independent of for large and fixed rj. [In practice, Monte Carlo simulations show 
that f{ri) is independent of A^ for A^ > 100]. 
We get in addition, 

<U>= ~3NT [1 - fir])] . 
In the dilute limit, rj ^ and we find the perfect gas value 

/(O) = 1 . 
Equating eqs.( [4.12[ ) and ( [4.14| ) yields, 

rv^l- fix) 



<!>Niv) = r dx 
Jo 



X 



Relevant thermodynamic quantities can be expressed in terms of the function fit]). We 
find for the free energy from eq. 



F = ~3NT\og\J—]-3NT P dx^-^^ . (4.15) 
\ V 27r / JO X 

We find for the total energy 

E = -3NT[^ - fir^)] , 

for the chemical potential. 



and for the entropy 



V 



_|iV + 3Ariog(J^| +3Ar r dx^-^^ + 3N firj) . (4.16) 
2 \ y In Jo X 



The specific heat at constant volume takes the form |l[| 

T fdS\ 



V 



fiv) 



V fiv) - \ 



(4.17) 



where we used eg. ( [4. 161) . This quantity is also related to the fluctuations of the potential 
energy (A?7)^ and it is positive deflned, 



cv 



Here, 



{Auy 



<U^ > - <U 



3[/(r/)-r]/'(r])-l] 



(4.18) 



The speciflc heat at constant pressure is given by IT|] 



Cp = Cv 



T 




2 
V 


N 







and then, 



Cp = Cv + 

3 

= -2' 



[fiv)-vnv)? 

fiv) + hf'iv) 

4/(r/) ifiv)-vf'iv)) 



fiv) + Ivf'iv) 

The isotherm (-R't) and adiabatic iKs) compressibilities take the form 



(4.19) 



Kt = — 



V 



dp, 



V 



NTfir,) + hn7^) 



K, = —. 



V 



dp, 



Cv 



Cp 



The speed of sound Vs can be written here as ||T6[ 



Cp V 

' cvN 



' dp ' 
dV 



N 



T 



Ncv 



' dp ' 



dp ' 



Therefore, 



T 



Ifiv) - vfiv)f 



fiv) - vf'iv) 



+ fiv) 



^vf'iv) ■ 



(4.20) 



We see that the large N hmit of the self-gravitating gas is special. Energy, free energy 
and entropy are extensive in the sense that they are proportional to the number of particles 
N (for flxed rj). They all depend on the variable t] = ^^^^^ which is to be kept flxed for 
the thermodynamic limit {N — > 00 and V — > 00) to exist. Notice that rj contains the ratio 
N/L = N V~^^^ and it is not an intensive variable in the usual sense. Here, the presence of 
long-range gravitational situations calls for a new type of variables in the thermodynamic 
limit. 
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A. Short-distance cutoff 

At short distance the particle interaction for the self-gravitating gas in physical situations 
is not gravitational. Its exact nature depends on the problem under consideration (opacity 
limit, Van der Waals forces for molecules etc.). We shall just assume a repulsive short 
distance potential. That is, 

for r > a 

for r < a (4.21) 

where r = \ri — fj\ stands for the distance between the particles and a << 1 is the short 
distance cut-off. 

The presence of the repulsive short-distance interaction prevents the collapse (here un- 
physical) of the self-gravitating gas. In the situations we are interested to describe (inter- 
stellar medium, galaxy distributions) the collapse situation is unphysical. 



Va{r) = — 
r 

Va{r) = +- 



B. The diluted regime: r/ « 1 

We can obtain the thermodynamic quantities as a series in powers of r] just expanding 
the exponent in the integrand of ^n{v) [eg - ([4.91) ]. 
To first order in rj we get, 

rl A 



rl fl 

Jo Jo f-Ji 
2 Jo Jo \ri — 

3iN - 1) bo V + Oia^) + 0(7]^) . (4.22) 



where the coefficient Bq is just a pure number. For the cubic geometry chosen, it takes the 
value 

bo = - / {l-x)dx il-y)dy \ = 0.31372... . 

o JO Jo Jo va;^ + 1/^ + 2;^ 

To first order in r] we see that the cutoff effect is negligeable ~ 0{a^). 

We therefore find in the low density limit using eqs.( [4.12D , ([4.14|) and ( [4.22| ) 

^ = f{^) = l-bov + 0{v'), (4.23) 

in the large limit. 

Furthermore, the speed of sound approaches hnearly in rj to its perfect gas value. 
We used here eqs.(^:20D and ( ^1231) . 
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C. Monte Carlo calculations 



We have applied the standard Metropohs algorithm to the self-gravitating gas in a cube 
of size L at temperature T. We computed in this way the pressure, the potential energy 
fluctuations and the average particle distance as functions of rj. 

Two different phases show up: for t] < rjc "we have a non-perfect gas and for t] > rjc it 
is a condensed system with negative pressure. The transition between the two phases is 
very sharp. We can say that the phase transition is of first order since there is a jump in 
the entropy. However, it is an unusual phase transition since the pressure is negative in the 
denser phase. In particular, the phases cannot coexist since the pressure has opposite sign 
in the two phases. This phase transition is associated with the Jeans unstability. 

We plot in figs. 1-3 /(r/) = pV/[NT] , (Af/)^ and the speed of sound squared v^/T as 
functions of f]. 

We find that for small rj, the Monte Carlo results for pV/[NT] well reproduce the ana- 
lytical formula ( [4.23| ). pV/[NT] monotonically decreases with rj. When rj gets close to the 
value r^c ~ 1-6 a phase transition suddenly hapens andpl^/[A^T] becomes large and negative. 

< r > monotonically decreases with t] too. Near rjc ,< r > has a sharp decrease. 

In the Monte Carlo simulations the phase transition to the condensed phase happens for 
rj = TjT slightly below 7]^. For t]t < f] < f]c, the gaseous phase may only exist as a metastable 
state. We find that rjc — rjT decreases with the number of particles in the simulation. For 
example, fjc — Vt ^ 0.2 for = 2000. 

Both, the values of pV/[NT] and < r > in the condensed phase depend on the cutoff a. 
The Monte Carlo results for rj > rjc can be approximated as 

pV Tl 

-— = f{r])^l--J- , <r>~1.5a. 
Nl K a 

where K ~ 30. Therefore, the latent heat of the transition is 

7] 



1 - 



< 



K a\ 

The behaviour of pV/[NT] near rjc in the gaseous phase can be well reproduced by 



NT 



fiv) "^^ A (r^, - (4.24) 



where A ~ 0.68, ri^ ~ 1.59 and B ~ 0.22. 

In addition, the behaviour of (A[/)^ in the same region is well reproduced by 

{AUy "I?^ C{r],-r]f-' (4.25) 

with C ~ 0.2 and with the exponent 5 — 1 as it must be since (Af/)^ grows as /'(r/) [see 
eq. ([4.18|) ]. [Notice that for finite A^, (Af/)^ will be finite albeit very large at the phase 
transition]. 

We thus find a critical region just below t]c where the energy fluctuations tend to infinity 
as ?7 t Vc- 
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Vg, cy and Ks turn to be positive in the whole interval < r] < rjc while cp and Kt are 
positive for < ?7 < ?7o and change sign (and diverge) at the point r^o < rj^ defined by 

fiVo) + lvo fivo) = . 

We find rjQ ~ 1.49 from the Monte Carlo simulations. We find that cp and are negative 
for riQ < T] < Tjc- 

The specific heat behaviour near the transition follows from eqs.( |4.17 ), ( 4.19| ) and ( |4.24 ), 



Cp 



+ (ric-vY 



cv 



3ABr]c{r]c- rf) 



B-l 



+ 



B 



That is, while cp tends to a negative value, cy grows without bound when r] | 7]^. The 
usual thermodynamic inequalities |]T| forbiding such negative values do not apply here due 
to the inhomogeneity of the self-gravitating gas at thermal equilibrium. 

We find for the speed of sound near the phase transition from eqs. ([4.20|) and ( [4.24|) , 



T 



6 



+ 



That is, the speed of sound tends to a constant value when r] | rjc- 

In the condensed phase, becomes negative indicating that there is no sound propaga- 
tion in such state. 

We verified that the Monte Carlo results in the gaseous phase {rj < rjc) are cutoff inde- 
pendent for 0.001 >a> 0.0. 



D. Particle Distribution 

The particle distribution at thermal equilibrium obtained through the Monte Carlo sim- 
ulations is inhomogenous both in the gaseous and condensed phases. 

In the gaseous phase we find from the Monte Carlo values of the particle density distri- 
bution that the mass within a volume V = scales as 

M = CR^ (4.26) 

where C is a i? independent constant and D takes values in the range, 

D = 1.9- 2.2 . 

near the phase transition. That is for lA < rj < rjc. For smaller rj, D increases towards the 
value D = 3 for rj ^ 0. That is, the distribution becomes uniform in the perfect gas limit, 
as expected. 

The exponent D found here suggest the presence of a fractal distribution near the the 
critical point pV/[NT] = O"*" for a very large number N of particles. 
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V. DISCUSSION 



We present here a set of new results for the self-gravitating thermal gas obtained by Monte 
Carlo and analytic methods. In particular, they confirm the general picture of the thermal 
self-gravitating gas. Namely, a gaseous phase for higher temperature and lower density and 
a condensed phase for lower temperature and higher density |[10|JT3|JT^ . Actually, we find 



more appropriate to characterize the phases by the sign of the pressure: positive pressure in 
the gaseous phase and negative pressure in the condensed phase. The pressure plays here 
the role of order parameter. 

The parameter t] [introduced in eg. ([4. 61)] can be related to the Jeans length of the system 




dj = \ 77^=' (5-1) 



where p = N/V stands for the number volume density. Combining eqs.( [4.6| ) and ( |5.1| ) yields 



We see that the phase transition takes place for dj ~ L. [The precise numerical value of the 
proportionality coefficient depends on the geometry]. For dj > L we find the gaseous phase 
and for dj < L the system condenses as expected. 



Contrary to mean field treatments [12,10|, we do not assume here an equation of state 



but we obtain the equation of state for the canonical ensemble [see eq.( |4.14| )]. We find at the 
same time that the relevant variable is here t] = Gm'^N/[V^^^T]. The relevance of the ratio 
Gm? /[V^^^T] has been noticed on dimensionality grounds |T^. However, dimensionality 
arguments alone cannot single out the crucial factor N in the variable t]. 

The crucial point is that the thermodynamic limit exist if we let ^ oo and V oo 
keeping i] fixed. Notice that i] contains the ratio V~^^^ and not N/V. This means that 
in this thermodynamic limit V grows as A^^ and thus the volume density p = N/V decreases 
as ~ A^~^. f] is to be kept fixed for a thermodynamic limit to exist in the same way as the 
temperature. pV, the energy E, the free energy, the entropy are functions of t] and T times 
A^. The chemical potential, specific heat, etc. are just functions of rj and T. 

The divergent growth of the energy fiuctuations (A[/)^ near the phase transition has 
been previously noticed []TU|,|T3|,|T1[ . We find here the precise behaviour of (Af/)^ for 77 t r/c 
using Monte Carlo methods [eq. ( |4.25| )]. 

In refs. [@-§l we worked in the grand canonical ensemble around the point where 
log Zgc = 0. This precisely corresponds to pV/[NT] = which is the critical point in 
the canonical treatment given here. The presence of a critical region where scaling holds 
supports the previous work in the grand canonical ensemble 
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FIG. 3. The speed of sound squared divided by the temperature, v^/T, as a function of 77 
in the gaseous phase from Monte Carlo simulations. Notice that v1/T takes the value 1/6 at the 
critical point r) = ijc- 
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